Multifractal fluctuations in the survival probability of an open quantum system 
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We predict a multifractal behaviour of transport in the deep quantum regime for the opened 
5— kicked rotor model. Our analysis focuses on intermediate and large scale correlations in the 
transport signal and generalizes previously found parametric mono-fractal fluctuations in the quan- 
tum survival probability on small scales. 
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I. INTRODUCTION 

Multifractal analysis of fluctuating signals is a widely 
applied method to characterize complexity on many 
scales in classical dynamics 0, or in the analysis of a 
given time series (without any a priori knowledge on the 
underlying dynamical system which generated the series) 


On the quantum level, multifractal behaviour was 
found in the scaling of eigenfunctions in solid-state trans- 
port problems [3j- As far as we know, there have been, 
however, very few attempts to use the method of multi- 
fractal analysis to directly characterize transport proper- 
ties such as conductance (across a solid state sample) or 
the survival probability (in open, decaying systems). Of- 
ten it is indirectly argued that the multifractal structure 
of the wave functions at critical points (at the crossover 
between the localized and the extended regime) imprints 
itself on the scaling of transport coefficients Other 
works found a fractal scaling of local transport quantities, 
such as hopping amplitudes 5] or two-point correlations 
|(|. At criticality [|| predicts, e.g., a multifractal scaling 
of the two-point conductance between two small interior 
probes within the transporting sample. 

In this paper, we directly study the fluctuations prop- 
erties of a global conductance like quantity in a regime 
of strong localization (Anderson or dynamical localiza- 
tion in our context of quantum dynamical systems). The 
studied quantity is the survival probability of an open, 
classically chaotic system, which in the deep quantum 
realm was found to obey a monofractal scaling if certain 
conditions on the quantum eigenvalue spectrum are ful- 
filled 0,13- In particular, the distribution of decay rates 
of the weakly opened system needs to obey a power-law 
with an exponent 7 ~ — 1, which translates into an an- 
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alytic prediction for the corresponding box counting di- 
mension (of the survival probability as a function of a 
proper scan parameter): Dbc — 2 — \j\/2. 

A more detailed, yet preliminary numerical analysis 
of the decay rate distribution for our model system (to 
be introduced below) has found that two scaling regions 
can be identified . While for small rates the probabil- 
ity density function p(T) scales as T -1 , at larger scales it 
turns to r~ 3 / 2 - which is expected for strongly transmit- 
ting channels from various models for transport through 
disordered systems Here we ask ourselves whether 
this prediction of a smooth variation in the scaling of the 
monofractal behaviour (induced by the smoothly chang- 
ing exponent 7) can be generalized to characterize the 
fluctuations on many scales using from the very begin- 
ning the technique of multifractal analysis. Before we 
present our findings on the multifractal scaling of the 
parametric fluctuations of the survival probability, we 
introduce the kicked rotor system and our numerical al- 
gorithm for the multifractal analysis in the subsequent 
two sections. 



II. OUR TRANSPORT MODEL AND THE 
CENTRAL OBSERVABLE 

The 8— kicked rotor is a widely studied, paradigmatic 
toy model of classical and quantum dynamical theory 
[T2I IT^ . Using either cold or ultracold atomic gases, 
the kicked rotor is realised experimentally by preparing 
a cloud of atoms with a small spread of initial momenta, 
which is then subjected to a one-dimensional optical lat- 
tice potential, flashed periodically in time In good 
approximation, the Hamiltonian for the experimental re- 
alization of the rotor on the line (in one spatial dimen- 
sion) reads in dimensionless units |14|: 

2 00 
H(t') = — + k cosx ^<5(t' - tr) . (1) 
2 f=i 
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The derivation of the one-period quantum evolution op- 
erator exploits t he sp atial periodicity of the potential by 
Bloch's theorem [T^J. This defines quasimomentum (3 as 
a constant of the motion, the value of which is the frac- 
tional part of the physical momentum p in dimensionlcss 
units p = n + [3 (n 6 N). Since (3 is a conserved 
quantum number, p can be labelled using its integer part 
n only. The spatial coordinate is then substituted by 
6 = x mod(27r) and the quantum momentum operator 
by J\f = —id 1 89 with periodic boundary conditions. The 
one-kick quantum propagation operator for a fixed f3 is 
thus given by 



Ur = e 



ik cos(e) e -ir(Ar+0) 2 / 2 



(2) 



In close analogy to the transport problem across a 
solid-state sample, we follow @ OH to define the quan- 
tum survival probability as the fraction of the atomic 
ensemble which stays within a specified region of mo- 
menta while applying absorbing boundary conditions at 
the "sample" edges. If we call ip(n) the wave function in 
momentum space and ri\ < ni the edges of the system, 
absorbing boundary conditions are implemented by set- 
ting ij)(n) = if n < ri\ = —1 or n > ri2 = 251. This 
truncation is carried out after each kick, and it mimics 
the escape of atoms out of the spatial region where the 
dynamics induced by the Hamiltonian Q takes place. If 
we denote by P the projection operator on the interval 
]ni, ni\ the survival probability after t kicks is: 



p SUI (t)= (p^)V(«,o) 



(3) 



The early studies of the fluctuation properties of P sur 
focused on its parametric dependence on the quasimo- 
mentum (3 [jj, UJi ■ While (3 is hard to control experimen- 
tally on a range of many scales (with a typical uncer- 
tainty of 0.1 in experiments with an initial ensemble of 
ultra-co\d atoms |lfijV some of us recently proposed to 
investigate the parametric fluctuations as a function of 
the kicking period t (see Eq. Q), which can be easily 
controlled on many scales in the experimental realization 
of the model even with laser-cooled (just "cold") atoms 
0. 

In Figure 1 we present the survival probability P sur 
of the opened kicked rotor in the deep quantum regime 
(i.e., at kicking periods t = h c g > 1 [Tj]) as a func- 
tion of the two different scan parameters (3 and r. The 
global oscillation with a period of the order 1 in Fig. 1 (a) 
originates from the /3-dependent phase term J\f(3 in the 
evolution operator J2J, and can be understood qualita- 
tively by remembering the Bloch band structure of the 
corresponding quasienergy spectrum as a function of [3 
[lflj . No such oscillating trend is found for the graph as 
a function of the kicking period. Nonetheless, in the fol- 
lowing, we use a well developed variation of the standard 
multifractal algorithm, which intrinsically takes account 
of such global, yet irrelevant trends in the signal function 



Fluctuation Analysis (MF-DFA) 18| are now explained 
before we present our central results which indicate the 
multifractal scaling of data sets as the ones shown in Fig. 
1. 



III. MULTIFRACTAL DETRENDED 
FLUCTUATION ANALYSIS 



The MF-DFA is a gene ralization of the DFA method 
ori gina lly proposed by |17| , and it is extensively described 



m [Tj|. In the recent years it was used, for instance, 
to investigate the nonlinear properties of nonstationary 
series of wind speed recor ds fl^| , electro-cardiograms [2(| , 
and financial time series |2l| . 

The method consists of five steps. First the series 
{ x i}iLi is integrated to give the profile function: 



y 



(4) 



where x is the average value of Xi. The profile can be 
considered as a random walk, which makes a jump to the 
right if Xi — x, is positive or to the left side if X i X IS 
negative. In order to analyze the fluctuations, the profile 
is divided into N s = int(N/s) non-overlapping segments 
of length s, and, since usually N is not an integer multiple 
of s, to avoid the cutting of the last part of the series, 
the procedure is repeated backwards starting from the 
end to the beginning of the data set. In each segment v 
we subtract the local polynomial trend of order k and we 
compute the variance: 

f\ v , = {v[iy - 1 > + i] - y k M 2 , (5) 
s i=i 

for v = l,...,N s , and 



with v — N s + 1, . . . , 2A^ S for the backward direction. 
The order of the polynomial defines the order of the MF- 
DFA too, therefore we may speak about MF-DFA(l), 
MF-DFA(2), ... , MF-DFA(fc). 

The fourth step consists on the averaging of all seg- 
ments to obtain the q-th order fluctuation function for 
segments of size s: 



F»(s) 



1 



2N S -\ 1 /l 



(7) 



In the last step we determine the scaling behaviour of 
the fluctuation function by analyzing the log-log plots of 
F q (s) versus s for each value of q. If the series is long- 
range correlated F q (s) increases for large s as a power 
law: 



F q (s) 



Mq) 



(8) 
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Since the number of segments becomes too small for very 
large scales (s > N a /4), we usually exclude these scales 
for the fitting procedure to determine h(q). The MF-DFA 
reduces to the standard DFA for q = 2, while the scaling 
exponent h(q) can be related to the standard multifractal 
analysis considering stationary time series, in which h(2) 
is identical to the Hurst exponent H , therefore, h(q) can 
be considered a generalized Hurst exponent. Monofrac- 
tal series indeed show a very weak or no dependence of 
h(q) on q. By example, for monofractal series as white 
noise, the generalized Hurst exponent is H = 1/2 for all 
q. On the contrary, for multifractal time series, h{q) is 
a function of q and this dependence influences the mul- 
tifractality of the process. Referring to the formalism of 
the partition function: 



Z q (s) 



E 



Vvs - V{ v -V)s\ 



(9) 



where r(q) is the Renyi exponent, to which the h(q) is 
related by: 



T (q) = 1 - qh(q). 



(10) 



Now we are able to use the formalism of the multifractal 
spectrum [23 f(a) to characterize the data set: 



dr(q) 



h(q) 



dh{q) 



dq ^ * dq (11) 

/(a) = qa - r(q) = q[a - h(q)] + 1 . 

The generalized dimensions are expressed as a function 
of t(<?) or h(q) [H: 



D„ 



r(q) qh(q) - 1 



(12) 







which cannot be not straightforwardly defined for q 
and q = 1. 

If the signal is multifractal, the spectrum /(a) has ap- 
proximately the form of an inverted parabola. As signi- 
ficative parameters for its characterization we considered 
the point olm corresponding to the maximum of f(as), 
and its width W a considered for a fixed q interval. In 
other words, «m represents the a value at which is situ- 
ated the "statistically most significant part" of the time 
series (i.e, the subsets with maximum fractal dimension 
among all subsets of the series). The width W a is re- 
lated to the dependence on h(q) from q. The stronger 
this dependence, the wider is the fractal spectrum (cf., 
eq. mj). 



IV. RESULTS 

We performed a MF-DFA of order k = 1 on data sets 
produced by scanning the (3 or r parameter, respectively, 
over 10 5 data points, and considering different interac- 
tion times from t = 250 to t — 10000 kicks. The analysis 



performed with higher order (fc = 2 and 3) polynomial 
detrending for some of the series produced basically the 
same results. Furthermore, we tested our numerical al- 
gorithm on a monofractal time series (white noise) and 
a well known multifractal process (binomial multifractal 
model 24]). For these two test series we reproduce the 
known analytical resuls, with a precision better than 1%. 

A full analysis for t = 6000 (see Figure 1) is shown in 
Figures 2 and 3 for the (3 and r scanned series, respec- 
tively. Tables HI and ITU collect the multifractal parameters 
ctM and W a , which were computed for t — 250 . . . 10000. 
Analogously to [2{|, we defined W a as the width of the 
parabolic form of f(ct) between the points corresponding 
to q = — 3 and q = 3. 

Figure 2(a) shows the scaling behaviour of the fluctu- 
ation function F q (s), with q 6 [—5,5]. Here s represents 
the index of the scanning parameter (3, while the fit was 
performed in the zone log(s) 6 [1.6,2.7] (corresponding 
to s S [40,500]). In Figure 2(b) we report the depen- 
dence of h(q) on g, revealing the multifractal nature of 
the data set. In order to better characterize the mul- 
tifractality and to highlight how it changes among the 
different analyzed series, we have computed the MF spec- 
trum f(a) (cf. 2(c)). Figure 2(d) shows the variation of 
the multifractal parameters for the different interaction 
times considered. After a fast decrease, both the param- 
eters tend to converge around the values olm = 1-29 and 
W a = 0.2 (see also Tab. QJ. Very similar results were 
obtained for the r scanned series (cf., Figure 3 and Table 
ITT|) . Comparing the values of Tables UJ and [H] we can say 
that both the r and the (3 scanned series have essentially 
the same multifractality. 

Even if we cannot a priori predict the asymptotic sim- 
ilarity between the two series of r and /?, we can a pos- 
teriori interpret this result: both parameters enter not 
equally yet similarly in the phase of the second factor on 
the right of eq. (2). As a consequence, the restriction of 
[3 to the unit interval does make no difference to the, in 
principle, unboundedness of r (in fact, to avoid different 
dynamical properties of the system, r was chosen in a 
restricted window too, cf. 0)- 

In general, two types of multifractality can be distin- 
guished, and both of them require different scaling ex- 
ponents for small and large fluctuations. (I) The mul- 
tifractality can be due to the broad probability density 
function for the values, and (II) it can also be due to 
different long range correlations for small and large fluc- 
tuations. The simplest way to distinguish between the 
mentioned two cases is to perform the analysis on a ran- 
domly reshuffled series. The shuffling destroys all the 
correlations, and the series with multifractals of type (II) 
will exhibit a monofractal behaviour with h Bhut (q) = 0.5 
and W a — 0. On the contrary, multifractality of type (I) 
is not affected by the shuffling procedure. If both (I) and 
(II) are present the series will show a weaker multifrac- 
tality than the original one. 

We applied the shuffling procedure to the series showed 
in Figure 1. The procedure destroyed the multifractality 
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of both series since for both the sequences we obtained 
h(q) = 0.51 ± 0.01 for q e [-5,5]. The dependence on 
q was so weak that we were not able to compute any 
reliable f(a) spectrum, which, in this case, can be con- 
sidered singular, i.e., with W a ~ 0. 

V. CONCLUSIONS 

We studied the quantum kicked rotor, a paradigmatic 
model of quantum chaos, which describes the time evo- 
lution of cold atoms in periodically flashed optical lat- 
tices. Imposing absorbing boundary conditions allows 
one to probe the transport properties of the system, here 
expressed by the survival probability on a finite region 
in momentum space. For a fixed interaction time, the 
quantum survival probability depends sensitively on the 
parameters of the system, and our application of the de- 
trended multifractal method shows that clear signatures 
of a multifractal scaling of the survival probability are 
found, as either the kicking period or quasimomentum 
is scanned. Our results generalize the previously pre- 
dicted mono- fractal structure of the signal 0, IS EI , by 
characterizing long-range correlations in the parametric 
fluctuations. In agreement with the monotonic increase 
of the box counting dimension with the interaction time t 
and its saturation after t > 5000 observed in Q , we found 



a systematically decreasing value for the maximum olm 
of the MF spectrum and of its widths Wm ■ Both of these 
two values also tend to saturate for t > 5000. 

Future work along the lines of @ will be devoted to 
check in detail whether traces of the here predicted mul- 
tifractality could be observed under real-life experimen- 
tal conditions (e.g., for short interaction times and finite 
resolutions in the scanning parameter Q). 
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FIG. 1: The series P sttr (/3) (a) and P sur (r) (b) after an 
interaction time of t = 6000 kicks. Other paramters are k = 5, 
and r = 1.4 in (a) and f3 = in (b), respectively. Both 
sequences extend over 10 5 sampling points along the shown 
intervals. 




Kicks 



FIG. 2: (a) The decimal logarithm of the fluctuation function 
F q (s) for q G [—5, 5] for the f) scanned series after an interac- 
tion time of 6000 kicks. The fitting procedure was performed 
in the zone log(s) € [1.6, 2.7] (corresponding to s € [40, 500]). 
The curves were vertically shifted for better reading, (b) The 
spectrum of the generalized Hurst exponents h(q); its strong 
dependence on q indicates the multifractal behaviour. The 
error bars show the uncertainty arising from the fits to the 
curves in (a), (c) The f(a) spectrum with olm = 1-29 and 
W a = 0.39. The dotted lines indicate the a interval used to 
compute Wa- (d) The multifractal parameters ckm and W a 
as a function of the interaction time. After a strong, initial 
variation, the value «m shows a saturation towards the value 
cum ~ 1.3. The width W a shows approximately the same be- 
haviour, and tends to saturate towards the value W a w 0.2. 




FIG. 3: (a) The decimal logarithm of F q (s) with q £ [-5,5] 
for the r scanned series with t = 6000. The fitting procedure 
was performed in the zone log(s) £ [1.6, 2.7] (corresponding 
to s 6 [40,500]). (b) The spectrum of the generalized Hurst 
exponents h(q). (c) The f{a) spectrum with ckm = 1-29 and 
W a — 0.28. (d) The variation of the multifractal parameters 
«m and W a with the interaction time, om shows a saturation 
towards the value au ~ 1-3, while W a tends to saturate 
around the value W a ^0.2. 
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TABLE I: Multifractal parameters of the /3-scanned data sets 
for different interaction times. The estimated error due to the 
fitting procedure described in section ITTT1 is about ±0.02 for 
cum and ±0.05 for W a - 

OtM W a 



t=250 1.63 1.19 

t=500 1.45 0.85 

t=1000 1.34 0.58 

t=2000 1.32 0.37 

t=4000 1.30 0.46 

t=6000 1.29 0.39 

t=8000 1.29 0.26 

t=10000 1.29 0.20 



TABLE II: Multifractal parameters of the r-scanned data sets 
for different interaction times (error estimates as stated for 
the preceding table). 

olm W a 



t=250 1.70 1.13 

t=500 1.48 1.12 

t=1000 1.35 0.9 

t=2000 1.33 0.52 

t=4000 1.30 0.44 

t=6000 1.29 0.28 

t=8000 1.27 0.19 

t=10000 1.27 0.20 



